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Based on angle resolved photoemission spectra measured on different systems at different dopings, 
momenta and photon energies, we show that the anomalously large spectral linewidth in the (it, 0) 
region of optimal doped and underdoped r^S^CaC^Os+s has significant contributions from the 
bilayer splitting, and that the scattering rate in this region is considerably smaller than previously 
estimated. This new picture of the electronic excitation near (n, 0) puts additional experimental 
constraints on various microscopic theories and data analysis. 

PACS numbers: 71.18.+y, 74.72.Hs, 79.60.Bm 



Angle resolved photoemission spectroscopy (ARPES) 
data from the (ir, 0) region of P^S^CaC^Og+a (Bi2212) 
have been one of the most important sources of infor- 
mation about the electronic structure of the high tem- 
perature superconductors (HTSC) jjj. The normal state 
spectra are very broad, with widths much larger than 
those from the nodal region (near (ir/2, tt/2)), indicating 
a large anisotropy in the scattering rate along the Fermi 
surface [g,|| . This anisotropy has been considered in var- 
ious theories that describe the anomalous transport and 
optical properties in the cuprates |4[|. In addition, the in- 
formation gathered from these spectra have helped to put 
additional important parameters into microscopic models 
0. On the other hand, the superconducting state spec- 
tra contains the well known peak-dip- hump structure m . 
The position of the dip was suggested to be related to the 
neutron (ir,ir) resonance mode [[7J, resulting in modelling 
of the tunnelling and ARPES data ||. The peak intensity 
in the peak-dip-hump structure has been interpreted as 
being related to the condensate fraction M , as discussed 
by various theories pOUO]. These studies constitute a 
significant part of the HTSC literature. 

Bi2212 has two coupled Cu02 planes in the unit 
cell and therefore bilayer splitting is naturally expected. 
However, it has been largely ignored in the studies men- 
tioned above, partly because of earlier reports of its ab- 
sence in the ARPES spectra J12],[l3| ■ Recently, this long- 
sought bilayer splitting was finally observed in overdoped 
Bi2212 [|14|-|l6fl . The two originally degenerate bands (one 
for each Cu02 layer) are split into bonding and antibond- 
ing bands due to the intra-bilayer coupling. In the (ir, 0) 
region, the amplitude of the bilayer splitting is found to 
be about 100 meV, comparable to the size of the super- 
conducting gap and the normal state band dispersion. 



As a result, the bilayer splitting causes a peak-dip-hump 
structure even in the normal state of heavily overdoped 
Bi2212 Jl^JT^pt, demonstrating that the intra-bilayer 
coupling plays an important role in the electronic struc- 
ture of the overdoped regime and should be seriously con- 
sidered in relevant theories. These results naturally raise 
the question of whether bilayer splitting exists in the op- 
timal and under-doped regimes where most experiments 
and analyses were conducted, and if it does, how it af- 
fects our understanding of the nature of the electronic 
excitations near (7r,0). 

In this letter, we report ARPES spectra from Bi2212 
and Bi2Sr2Cu06+a (Bi2201) for various dopings and pho- 
ton energies (hi/). The lineshapes of Bi2201 and Bi2212 
are similar in the nodal region, but very different near 
(ir, 0). In addition, Bi2212 spectra from the (tv, 0) region 
are strongly modified by hi/, in contrast to the weak pho- 
ton energy dependence of the Bi2201 spectra. We show 
that these results can be well explained by the underly- 
ing bilayer splitting effects in under and optimally doped 
Bi2212 and that the broad linewidth near (tt, 0) is, in 
large part, due to the bilayer splitting. These results are 
very different from the current, commonly- accepted pic- 
ture of the electronic excitations near (ir, 0), and there- 
fore requires the reexamination of many existing theories, 
and puts strong constraints on future theoretical models 
and data analysis. 

High quality Bi2212 and Bi2201 single crystals were 
grown by the floating zone technique. Bi2212 samples 
are labeled by the superconducting phase transition tem- 
perature T c of the sample with the prefix UD for under- 
doped, OP for optimally doped, and OD for overdoped. 
Bi2201 samples are labeled in the same way but in low- 
ercase. For example, UD83 represents a T C =83K under- 
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FIG. 1. (a) Bilayer-split Fermi surfaces of heavily over- 
doped OD65; the two weaker features are their superstructure 
counter parts. Solid and dashed lines represent the bond- 
ing and antibonding Fermi surfaces, respectively, (b) Normal 
state photoemission spectra of Bi2212 taken at (n, 0) for three 
different doping levels. Data were taken with hv — 22.7eV 
photon. Bars indicate identified feature positions, and trian- 
gles indicate possible feature positions. 



doped Bi2212 sample, while odl7 represents a T C =17K 
overdoped Bi2201 sample. Samples with Pb doping are 
labeled with the prefix "Pb" , except od33, which is doped 
with both Pb and La. The superconducting transi- 
tion widths, AT C , were less than 3K for all the samples 
used. Angle resolved photoemission experiments were 
performed at a normal incidence monochromator (NIM) 
beamline of the Stanford Synchrotron Radiation Labo- 
ratory, where the intensity of the second order light is 
extremely weak. Data were taken with a Scienta SES200 
electron analyzer with the angular resolution of 0.3 x 0.5 
degrees unless specified otherwise. The overall energy 
resolution varied from 10 meV to 18 meV at different 
hf's. This variation of the energy resolution does not af- 
fect any of our conclusions since the energy scales of the 
discussed features are much larger. The chamber pres- 
sure was better than 5 x 10 _ ll torr, and sample aging 
effects were negligible during the measurements. Un- 
less otherwise specified, normal state data were taken 
10-20K above T c . 

The most obvious signature of the bilayer splitting is 
double features in the Fermi surface and energy distribu- 
tion curves (EDC), as shown in Fig. 1 for heavily over- 
doped Bi2212 |L4J. Due to the anisotropic nature of the 
intrabilayer coupling, the amplitude of the bilayer split- 
ting is also anisotropic with zero splitting in the nodal 
region and maximum splitting at (tt, 0) [Fig. 1(a)]. The 
normal state (tt, 0) EDC of heavily overdoped Pb-OD65 
[Fig. 1(b)] clearly shows two features that exhibit a nor- 
mal state peak-dip-hump structure, and are assigned to 
the two bilayer split bands |pT[-|l6|. This was not ob- 
served in previous measurements on overdoped samples, 
mainly due to extrinsic factors such as energy and an- 
gular resolution. With a slight decrease of the doping 
(Pb-OD72), the two components of the (tt, 0) spectrum 




OD63,T=75K 
0.6 -0.4 




OD63,T=75K 

Pb-odl7.T=20K (*•' 

■0.6 -0.4 -0.2 



Energy relative to E F (eV) 



FIG. 2. Normal state spectra taken at (a) the d-wave node 
region, and (b) the (tt, 0) region for both Bi2212 and Bi2201 
at various dopings. The solid and dashed lines in the inset 
of panel (a) indicate the bonding and antibonding Fermi sur- 
faces of OD63 respectively, and the black dots indicate the 
momentum locations of the spectra. Data were taken with 
hv — 22.7 'eV . Angular resolution was 0.3° for OP90 and 
0.12° for others in panel a). 



are barely distinguishable. Compared with Pb-OD65, the 
two features become broader and their intensities smaller. 
For OP90, the spectrum is intrinsically too broad to dis- 
tinguish the two split features, which makes the identifi- 
cation of the bilayer splitting very difficult in this man- 
ner. 

To clarify this further, we looked for other signs of bi- 
layer splitting by comparing the spectra of Bi2212 with 
those of single layer Bi2201 at similar doping levels. We 
chose two pairs of samples: OP90 and op33, and OD63 
and odl7. Based on the empirical T c vs. doping for- 
mula p9| , they have doping levels of 0.16, 0.16, 0.22, and 
0.24, respectively. For the spectra taken in the nodal re- 
gion shown in Fig. 2a, Bi2201 and Bi2212 have similar 
lineshapes, and the linewidth varies only slightly for dif- 
ferent systems and experimental conditions. This holds 
true even for the heavily overdoped Bi2201 sample with a 
T c < 1.7K (doping level ~ 0.28). The situation is very dif- 
ferent for the spectra taken in the (tt, 0) region (Fig. 2b). 
For OD63, the spectrum consists of both the bonding 
and antibonding bands, while the spectrum at (tt, 0.2tt) 
mostly consists of the bonding band, because the anti- 
bonding band is above E F pi. We find that the Bi2212 
and Bi2201 spectra match at (7t,0.27t) almost perfectly, 
while those at at (tt, 0) do not because of the presence 
of the bonding band at higher energies. As far as the 
neai-Ep features are concerned, the spectra from both 
odl7 and OD63 have very similar linewidths at similar 
binding energies and momenta. This similarity between 
the OD63/odl7 low energy spectra can be attributed to 
their similar doping levels in each Cu02 plane. The 
OP90/op33 (tt, 0) spectra show a large mismatch simi- 
lar to the OD63/odl7 case, which can be naturally at- 
tributed to the additional spectral weight from the bond- 
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FIG. 3. (a) Cartoon of the initial and final state symme- 
tries along the c-axis for the photoemission process with the 
presence of the bilayer splitting, (b) The calculated c-axis 
contributions to the photoemission matrix elements for both 
bonding and antibonding state as functions of hv, as described 
in the text. 



ing band of OP90. On the other hand, without bilayer 
splitting (or intra-bilayer coupling) in OP90, properties 
of the Cu02 planes of OP90 and op33 should be simi- 
lar. It is then difficult to explain why the linewidths of 
OP90 and op33 are so dramatically different in the (ir, 0) 
region, considering that Bi2201 and Bi2212 are very sim- 
ilar in many other aspects such as the phase diagram, 
Fermi surface shapes, dispersion energy scales p0| , pl[ , 
and particularly, residual resistivity, which indicates the 
scattering caused by defects and impurities. The larger 
linewidth of op33, compared to odl7, may be attributed 
to enhanced correlation effects with decreased doping, 
presumably (jr, n) scattering due to increased antiferro- 
magnetic fluctuations |g,[||. We note that these normal 
state spectra were taken at different temperatures. How- 
ever, the thermal broadening is negligible compared to 
the peak widths, within the experimental temperature 
range. 

The above comparison between spectra from Bi2201 
and Bi2212 suggests the possible presence of bilayer split- 
ting in optimally doped samples. This is further sup- 
ported by photon energy dependence studies. As de- 
picted in Fig. 3(a), the antibonding and bonding states 
have opposite symmetry along the c-axis with respect to 
the midpoint between the two CuOi planes. As a con- 
sequence, their photoemission matrix elements respond 
differently to various experimental parameters, including 
the photon energy. Upon tuning hv, the spectral weight 
from the bonding and antibonding states will vary differ- 
ently, thus changing the overall spectral lineshapes. This 
can be further illustrated by an analysis of the photoe- 
mission matrix elements. Although comprehensive calcu- 
lations of the photoemission matrix element are still not 
feasible because of the complexity in the crystal structure 
and the photoemission process, as well as the electron- 
electron correlations, with reasonable assumptions and 
simplifications one can still study its behaviors on a qual- 
itative level, which turns out to be very helpful for the 
interpretation of the data on various occasions [£2| . 

We consider two 2D systems coupled via a certain 



bilayer interaction. The photoemission intensity for 
such a system with non-interacting electrons is I oc 
S a = a b M£,A a (k, to) where k, uj and A a are the momen- 
tum, energy for the final state, and the spectral func- 
tion respectively, while a and b represent antibonding 
and bonding bands. In the one-electron matrix element 
M a = (ipf\A. ■ pIV'm), ipf and ipi a are the final state and 
initial state single electron wave functions, A is the vec- 
tor potential of the photon field, and p = —ihV. As- 
suming ipi a (x,y, z) — 4>i{x,y)x a {z) and a free electron 
final state ipf(x,y,z) = e ? k * x +' lk vy+ lk * z ) the matrix ele- 
ment can be separated into the in-plane contribution My 
and the out-of-plane contribution M za . Under the dipole 
approximation, 

M a oc M\\ + M za 

= A\\/A z {e ik ^ +ik yy\ n \U^y)) + (e** z MXa(*)> 

where k x and k y are fixed to be {-k, 0) for both the ini- 
tial and final states. The first term contributes equally 
to both bonding and antibonding states, and the ratio of 
polarization, A\\/A z , is approximately constant in the ex- 
periment. Therefore, we can focus on M za as a function 
of hv. To further simplify, we assume 
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where "-" and "+" signs are for a = a and b, re- 
spectively, ^o is the intrabilayer distance, and (3 is an 
adjustable parameter reflecting how the electron wave- 
function is localized within a CuOi layer and is as- 
sumed to be j3 = 1/6 in the calculation. For the fi- 
nal state, the free electron approximation gives k z = 
[2m*h~ 2 (hv — § + v a )-(kl + ky)]~2, where we choose the 
photoelectron effective mass m* to be the free electron 
mass, the work function $ = 4.3 eV, and inner potential 
t>o = 7 eV in the calculation |23|. M za calculated with 
these parameters and simplifications is shown for both 
antibonding and bonding states in Fig. 3(b). The M za 's 
for the bonding and antibonding state have almost oppo- 
site behaviors with hv, and changes quite dramatically in 
the studied hv range. This causes the overall lineshape of 
the Bi22I2 (ir, 0) spectrum to alter significantly with hv 
as the relative weight of bonding/antibonding states os- 
cillates. In the case of optimally doped and underdoped 
systems, the centroid of the broad feature will shift. 

This is indeed observed in OD65, where the bilayer 
splitting has been clearly identified |L4[| . Fig. 4(a-c) show 
ARPES intensity taken in the (ir, 0) region at different 
hv's as a function of momentum and binding energy. Be- 
cause the NIM gives extremely weak second-order light, it 
is possible to directly compare spectra taken at different 
hv's. One can see that the relative intensities of the anti- 
bonding band (A) and bonding band (B) change with hv. 
At some photon energies, only one feature is prominent, 
while in others, both features are clearly visible. EDC's 
of OD65 at (tt,0) are plotted in Fig. 4(d) PJ. While 
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FIG. 4. (a-c) Normal state photoemission intensity as a function of momentum and binding energy in grayscale maps near 
the (-7T, 0) region (k x — n) of OD65 taken with hv = 18.4 eV, 20 eV and 32 eV respectively. The thick white lines indicate the 
dispersion of the bonding (B) and antibonding (A) bands, and lines are shown in dashed form when the feature is weak, (d-h) 
Normal state EDO's taken at (tt,0) for samples (d) OD65, (e) OP90, (f) UD83, (g) UD73, and (h) op33. at various hv's. Bars 
serve as guides for the centroids of spectral features. Note that the features indicated by bars near Ef may be Fermi cutoffs 
instead of real features. The separation between bars is not necessarily the splitting energy. 



some k z dispersion may exist, the data show strong bi- 
layer matrix element effects. One clearly sees that the 
relative intensities of the bonding and antibonding fea- 
tures vary drastically with hv. For optimally doped [Fig. 
4(e)] and underdoped Bi2212 [Fig. 4(f-g)], one does not 
see two clearly separated features. However, one can see 
the strong variation of the lineshape, and changes in the 
centroid of the feature. Although there arc some de- 
tailed variations from sample to sample, the spectra of 
underdoped and optimally doped Bi2212 change with a 
similar trend as the OD65. On the other hand, for the 
optimally doped single layer system Bi2201 [Fig. 4(h)], 
the peak position and the overall lineshape show virtually 
no photon energy dependence. The high binding energy 
background of op33 is a smooth function of binding en- 
ergy and photon energy. These indicate that the strong 
photon energy dependence of the Bi2212 spectra is due 
to the bilayer splitting. 

Because of the previous lack of evidence for bilayer 
splitting in optimally doped and underdoped Bi2212 
p2| , p3| , many analyses and calculations assumed its ab- 
sence. For example, the momentum distribution curves 
in this region were usually fitted by one Lorentzian H^ , 
when in fact it consisted of two Lorentzians separated by 
the bilayer splitting in momentum space. Bi2212 (tt, 0) 
spectra were discussed J^Land particularly, fitted with 
a one-component formula £6|. We show that even with 
a two-component model, there are various uncertainties 
involved in fitting the spectra, because the bonding and 
antibonding features are weighted by unknown factors at 
certain hv's, and generally too broad to be reliably sep- 
arated. For overdoped Bi2212, the dip position will be 
shifted by the overlapping of the two humps and the the 
superconducting peaks. 

Intra-bilayer coupling was assumed in some theories to 
explain the different temperature dependence behavior 
of c-axis and in-plane transport and optical properties of 
bilayer systems jg7|. Our results reinforce the assump- 
tions of these theories. On the other hand, we find that 



the quasiparticle in the (7r, 0) region of optimally doped 
Bi2212 should be similar to that of Bi2201, and thus 
much better defined than previously believed from ear- 
lier Bi2212 data. The quasiparticle lifetime is more than 
100% longer than obtained from previous EDC analyses, 
i.e., the scattering rate in this (w, 0) "hot spot" is not as 
large as previously believed, although an anisotropy of 
the scattering rate still exists in the optimally doped and 
underdoped regime, as is observed in Bi2201 ||j||. 
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